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Abstract 

We present new and improved data structures that answer 
exact node-to-node distance queries in planar graphs. Such 
data structures are also known as distance oracles. For any 
directed planar graph on n nodes with non-negative lengths 
we obtain the following]]] 

• Given a desired space allocation S 6 [nlglgn,n 2 ], we 
show how to construct in O(S) time a data structure 

of size O(S) that answers distance queries in 0(n/\/~S) 
time per query. The best distance oracles for planar 
graphs until the current work are due to Cabello (SODA 
2006), Chen and Xu (STOC 2000), Djidjev (WG 1996), 
and Fakcharoenphol and Rao (FOCS 2001). For a £ 
(1,4/3) and space S = n CT , we essentially improve the 
query time from n 2 / S to \Jn 2 /S. 

• As a consequence, we obtain an improvement over the 
fastest algorithm for fc-many distances in planar graphs 
whenever k 6 [\/n, n) . 

• We provide a linear-space exact distance oracle for pla- 
nar graphs with query time 0(n 1 ^ 2+e ) for any constant 
e > 0. This is the first such data structure with prov- 
able sublinear query time. 

• For edge lengths > 1, we provide an exact distance ora- 
cle of space O(n) such that for any pair of nodes at dis- 
tance £ the query time is 0(min{£, ^/n}). Comparable 
query performance had been observed experimentally 
but could not be explained theoretically. 

Our data structures with superlinear space are based on the 
following new tool: given a non-self-crossing cycle C with 

c = 0{s/n) nodes, we can preprocess G in O(n) time to 
produce a data structure of size 0(n lglgc) that can answer 

the following queries in O(c) time: for a query node u, 
output the distance from u to all the nodes of C. This 
data structure builds on and provides an alternative for a 
related data structure of Klein (SODA 2005), which reports 
distances to the boundary of a face, rather than a cycle. 

1 Introduction. 

A fast shortest-path query data structure may be of use 
whenever an application needs to compute shortest-path 
distances between some but not necessarily all pairs of 
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nodes. Indeed, shortest-path query processing is an in- 
tegral part of many applications [SomlOj . in particular 
in Geographic Information Systems (GIS), intelligent 
transportation systems [JHR96] . These systems may 
help individuals in finding fast routes or they may also 
assist companies in improving fleet management, plant 
and facility layout, and supply chain management. A 
particular challenge for traffic information systems or 
public transportation systems is to process a vast num- 
ber of queries on-line while keeping the space require- 
ments as small as possible ZarfJS] . Low space consump- 
tion is obviously very important when a query algorithm 
is run on a system with heavily restricted memory such 
as a handheld device |GW05| but it is also important 
for systems with memory hierarchies [HMZ03, AT05 , 
where caching effects can have a significant impact on 
the query time. 

While many road networks are actually not exactly 
planar |EG08[ IAFGW10] , they still share many proper- 
ties with planar graphs; in particular, many road net- 
works appear to have small separators as well. For this 
reason, planar graphs are often used to model various 
transportation networks. 

In the following, we provide shortest-path query 
data structures (also known as distance oracles) for 
planar graphs for essentially any specified space re- 
quirement. Throughout the paper we assume the edge 
lengths to be non-negativej^] Our results extend re- 
sults of Cabello |Cab06j . Chen and Xu [CXOOj . and 
Fakcharoenphol and Rao |FR06j and extend one re- 
sult and improve upon another result of Djidjev |Dji96| . 
Consider Figure [l] for an illustration and comparison. 

Theorem 1.1. Let G be a directed planar graph on 
n vertices. For any value S in the range S G 
[n lglgn, n 2 ], there is a data structure with preprocess- 
ing time 0(S lg 3 nj lglgn) and space O(S) that answers 
distance queries in O {nS- 1 ' 2 lg 2 n lg 3/2 lg n) time per 
query. 



2 Our results apply to graphs with negative-length edges by 
using reduced lengths induced by a feasible price function |Joh77| . 
The current best bound for computing a feasible price-function in 
a planar graph is 0(nlg 2 ra/ lglgn) |MWN10I . 



As a corollary, we obtain the following result on k- 
many distances [FMS91] . improving upon an algorithm 
of Cabello |Cab06| . which runs in time 0((kn) 2 / 3 + 
n 4 / 3 ). Our result is an improvement for k = o(n). For 
the range roughly k 6 [s/n, n) , our algorithm is faster 
by a factor of 0((n/fc) 2 / 3 ). 

Theorem 1.2. Let G be a directed planar graph on n 
vertices. The distances between k = ^(n 1 / 2 lgn/lglgn) 
pairs of nodes (si, t±), (s2, ^2), ■ ■ ■ (sfcj £fc) can 6e com- 
pirfed m time 0((/cn) 2 / 3 (lgn) 7 / 3 (lglgn) 2 / 3 ). 

We also give a data structure that keeps the space 
requirements as small as possible, i.e. linear in the size 
of the input. This is the first linear-space data structure 
with provably sublinear query time for exact point- 
to-point shortest-path queries. Nussbaum jNusllj has 
simultaneously obtained a similar result. 

Theorem 1.3. For any directed planar graph G with 
non-negative arc lengths and for any constant e > 0, 
there is a data structure that supports exact distance 
queries in G with the following properties: the data 
structure can be created in time 0(n\gn), the space 
required is 0(n), and the query time is 0(n 1 / 2+e ). 

More generally, our data structure works for the 
range S € [n, nlglgn] (using a non-constant e, see de- 
tails in Section (31), where we slightly improve upon the 
tradeoff obtained by the construction of Fakcharoenphol 
and Rao [FR06 . Combined with Theorem |1.1[ we pro- 
vide efficient distance oracles for any space bound S £ 
[n, n 2 }. 

The main techniques we use are Frederickson's 
r-division |Fre87j . Fakcharoenphol and Rao's imple- 
mentation of Dijkstra's algorithm [FR06] . and Klein's 
Mult iple-So urce Shortest Paths (MSSP) data struc- 
ture |Kle05j , for which we propose a more space-efficient 
alternative (see Section [4] for a detailed comparison) as 
follows. 

Theorem 1.4. Given a directed planar graph G on n 
nodes and a simple cycle C with c — 0(^/n) nodes, there 
is an algorithm that preprocesses G in 0(n lg 3 n) time 
to produce a data structure of size 0(n lglgc) that can 
answer the following queries in 0{c lg 2 clglgc) time: for 
a query node u, output the distance from u to all the 
nodes of C. 

Since Klein's MSSP data structure has found nu- 
merous applications, we believe that our data structure 
could be a rather useful tool in other algorithms as wellj^j 



a There is at least one application for approximate distance 
oracles |Somlll . 



The shortest-path query problem has been investi- 
gated heavily from an experimental perspective as well. 
Experimental results suggest that query times propor- 
tional to the shortest-path length are possible in prac- 
tice using algorithms based on the so-called arc-flag 
technique |Lau041 IKMS051 IHKMSOQj g] Hilger, Kohler, 
Mdhring, and Schilling make a statement about the (ex- 
perimental) worst-case behavior of their method: 

In all cases, the search space of our arc- 
flag method is never larger than ten times 
the actual number of nodes on the shortest 
paths |HKMS09l Section 6]. 

We can now actually prove a similar statement: 

In (provably) all cases, the search space of our 
method is never larger than a polylogarithmic 
factor times the length i of the shortest paths. 

More precisely: 

Theorem 1.5. For any planar graph G with edge 
lengths > 1 there is an exact distance ora- 
cle using space O (nig nlglgn) with query time 
0(min \l lg 2 £\g\g£, \fn\g 2 n}) for any pair of nodes 
at distance I. The preprocessing time is bounded 
by 0(n 1+e ) for any constant e > 0. 

Note that our data structure has query time propor- 
tional to the path length. In fact, our algorithm main- 
tains a Bellman- Ford-type invariant: after iteration i, 
the distance represents the minimum path length among 
all paths on 0{i) edges — the correct distance is com- 
puted after time roughly proportional to the minimum 
number of edges on a shortest path but we can only 
guarantee correctness after time 0(t). If we may further 

assume that, for some constant e > 0, all s 1 paths 

of length at most (1 + e)dG(s,t) have f2(/ic?(s, £)) e dges 
(where hc(s,t) denotes the number of edges (hops) on 
a minimum-hop shortest-path) , then our data structure 
can be constructed to have query time proportional 
to the minimum number of edges on a shortest path 
0(h G (s,t)). This assumption essentially means that 
any s — t path with significantly fewer edges than the 
shortest path is much longer. 

Our main contributions can be summarized as fol- 
lows: (i) We improve the worst-case behavior of pre- 
viously known distance oracles with low space require- 
ments, providing a distance oracle with fast preprocess- 
ing algorithm that works for the whole tradeoff curve 

4 The preprocessing algorithm of this technique first partitions 
the graph into regions Vi and thereafter labels each edge e for all 
regions Vi with a boolean flag spi(e) indicating whether e lies on 
any shortest path to any node in Vi. At query time, only edges 
leading towards the target region need to be considered. 



(we also provide the first one with linear space and sub- 
linear query time), and (ii) we make an important step 
towards proving the behavior observed in practice. As 
our main tool, (Hi) we provide a more space-efficient 
multiple-source shortest-path data structure. 

1.1 Related Work. Shortest-path query processing 
for planar graphs has been studied extensively. In this 
section, we give a brief review of previous results. We 
focus on the space-query time tradeoff. See Figure [T] 
for a summary of known results in comparison with 
ours. The tradeoffs previously had not been illustrated 
by a space vs. query time plot as in Figure [TJ indeed, 
the illustration suggests that a data structure like the 
one described in our main theorem (Theorem |1.1[ ) was 
bound to exist. 
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Figure 1: Tradeoff of the Space [5] vs. the Query time 
[Q] for different shortest-path query data structures 
on a doubly logarithmic scale, ignoring constant and 
logarithmic factors. The upper line represents the 
Q = n 2 /S tradeoff (completely covered by oracles of 
Djidjev [Dji96| ; the oracles of Arikati, Chen, Chew, 
Das, Smid, and Zaroliagis [ACC + 96 cover the range 
5 G [n 3 / 2 ,n 2 ]; SSSP (5 = n) and APSP (5 = n 2 ) 
also lie on this line). The lower line represents the 
Q = n/\/~S tradeoff; the result of Djidjev |Dji96| covers 
the ran ge 5 € [n 4 / 3 ,n 3 / 2 ]; Chen and Xu [CXOOj and 
Cabello |Cab06j extend this to 5 € [n 4/3 ,n 2 ]. The data 
structure of Fakcharoenphol and Rao [FR06J covers the 
point S — n and Q = y/n. We extend their results to 
the full range. 

For exact shortest-path queries, the currently best 
result in terms of the tradeoff between space and query 
time is by Fakcharoenphol and Rao (FR06j . Their 



data structure of space 0(n) can be constructed in 
time O(n) and processes queries in time 0(\/n). The 
preprocessing time can be improved by a logarithmic 
factor [KMW10] ■ (Applications to distance oracles are 
not explicitly stated in |KMW10j .) 

Some distance oracles have better query times. 
Djidjev |Dji96 , improving upon earlier work of Fcuer- 
stein and Marchetti-Spaccamela |FMS91j . proves that, 
for any 5 € [n, n 2 } , there is an exact distance oracle with 
preprocessing time 0(5) (which increases to 0{n\fS) 
for 5 S [n,n 3 / 2 ]), space 0(5), and query time 0(n 2 /S). 
For a smaller range, he also proves that, for any 5 £ 
[n 4 ' 3 ,n 3 ' 2 ], there is an exact distance oracle with pre- 
processing time 0(nv5), space 0(5), and query time 
0(n/v5). Chen and Xu |CXOO] . extending the range, 
prove that, for any 5 £ [n 4 / 3 , rt 2 ], there is an exact dis- 
tance oracle using space 0(5) with preprocessing time 
0(nVS) and query time 0(n/VS). Cabello [CabOfij . 
mainly improving the preprocessing time, proves that, 
for any 5 € [n 4 / 3 ,n 2 ], there is an exact distance ora- 
cle with preprocessing time and space 0(5) and query 
time 0(n/\/S). Compared to Djidjev's construction, 
the query time is slower by a logarithmic factor but the 
range for 5 is larger and the preprocessing time is faster. 
Nussbaum [Nusllj provides a data structure that main- 
tains both the tradeoff from Dji96, CXOOJ and the fast 
preprocessing time from |Cab06| . Nussbaum also pro- 



vides a different data structure with a "clean" tradeoff of 
space 0(5) and query time 0(rt/v / 5), however spend- 
ing time 0(5 3//2 /^/n) in the preprocessing phase. In 
our construction, we sacrifice another root-logarithmic 
factor in the query time (compared to |Cab06j ) but we 
prove the bounds for essentially the whole range of 5; 
our preprocessing time remains 0(5). 

If constant query time is desired, storing a com- 
plete distance matrix is essentially the best solution 
we have. Wulff-Nilsen WNIOa recently improved the 
space requirements to o(n 2 ). If the space is restricted 
to linear, using the linear-time single-source shortest 
path algorithm of Henzinger, Klein, Rao, and Subrama- 
nian [HKRS97 is the fastest known for exact shortest- 
path queries until the current work (Theorem 1.3) 



Nussbaum [Nusllj has simultaneously obtained a result 
similar to Theorem 11.31 

Efficient data structures for shortest-path queries 
have also been devised for restricted classes of pla- 
nar graphs [DPZOOl ICX001 and for restricted types of 
queries |Epp99[ IKK06I ISch981 IKle05j . If approximate 
distances and shortest paths are sufficient, a better 
tradeoff with 0(n) space and O(l) query time has been 
achieved |Tho04l lKle02l lKle05l IKKSllj . 

Based on separators, geometric properties, and 
other characteristics such as highway structures, many 



efficient practical methods have been devised [GSSD08 , 
ISS05|. IBFSS07| . their time and space complexities 
are however difficult to analyze. Competitive worst- 
case bounds have been achieved under the assumption 
that actual road networks have small highway dimen- 
sion [AFGW101 lADF+11] . While our preprocessing al- 
gorithm (Theorem 1.5| runs in almost linear time, some 
of the problems that appear in the preprocessing stage 
of practical route planning methods have recently been 
proven to be NP-hard [EDDW091 lBCK+10j . 

2 Preliminaries. 

2.1 Recursive r division of Planar Graphs. Let 

G = (V, E) be a planar graph with |V| = n. Let Ep be 
a subset of E, and let P = (Vp, Ep) be the subgraph of 
G induced by Ep. P is called a piece of G. The nodes of 
Vp that are incident in G to nodes of V \ Vp are called 
the boundary nodes of P and denoted by dP. 

An r -division |Fre87] of G is a decomposition into 
0(n/r) edge-disjoint pieces, each with 0(r) nodes and 
0(^/r) boundary nodes. We use an r-division with the 
additional property that, in each piece, there exist a 
constant number of faces, called holes, such that every 
boundary node is incident to some hole. Such a decom- 
position can be found in C^nlgr + nr" 1 / 2 lgn) (WNlOb 
by applying Miller's cycle separator [Mil86j iteratively. 

We use this r-division recursively. Let the base of 
the recursion be level 0, and let the top of the recursion 
be level k. The entire graph G is defined to be the only 
piece at level k. The pieces of level i of the recursion 
are obtained by computing an r^-division for each level- 
(i + 1) piece. The notation suggests that we may use 
a different parameter r in the r-division at every level 
of the recursion. Indeed, using a non-uniform recursion 
is important in obtaining Theorem 1.3 For a level-i 
piece P, the level-(i — 1) pieces obtained by applying 
the i — division to P are called the subpieces of P. 

We stress that the classification of nodes of a piece 
at any level as boundary nodes is with respect to G (and 
not P). This implies that, if v is a boundary node of 
a level-i piece, then v is also a boundary node of any 
lower- level piece that contains v. This generalizes the 
decomposition used in FR06 , where Miller's separator 
is used at each level rather than an r-division. 

2.2 Klein's Multiple-Source Shortest Paths Al- 
gorithm. Klein [Kle05 gave a multiple-source shortest 
path (MSSP) algorithm with the following properties. 
The input consists of a directed planar embedded graph 
G with non-negative arc-lengths, and a face /. For each 
node u in turn on the boundary of /, the algorithm com- 
putes (an implicit representation of) the shortest-path 
tree rooted at u. This takes a total of O(nlgn) time 



and space. Subsequently, the distance between any pair 
(u, v) of nodes of G where u is on the boundary of / 
(and v is an arbitrary node), can be queried in O(lgn) 
time. In particular, given a set S of 0(*/n) nodes on the 
boundary of a single face, the algorithm can compute all 
S-to-S distances in 0{n\gn) time. 

We propose a more space-efficient alternative in Sec- 
tion |4j wherein we also provide a detailed comparison. 

2.3 Dense Distance Graphs and Efficient Im- 
plementation of Dijkstra's Algorithm. The dense 
distance graph for a piece P, denoted by DDGp, is 
the complete graph on dP, the boundary nodes of P, 
such that the length of an arc corresponds to the dis- 
tance (in P) between its endpoints. The dense dis- 
tance graph for all pieces P in the r-division can be 
computed in 0(|V(G)| lg | V(G)|) time and space using 
Klein's MSSP; For each piece P, compute MSSP data 
structures in 0(|P|lg|P|) time and space a constant 
number of times, specifying a different hole of P as the 
distinguished face at each run. Then query the MSSP 
data structures for the distances between the boundary 
nodes in 0(\dP\ 2 lg \P\) time. Since in an r-division 
\dP\ = y/\P\, this takes 0(\P\ lg |P|) time and space per 
piece, for a total of 0(\V(G)\\g \V(G)\) for all pieces. 

Let V be a set of pieces (not necessarily at the same 
level), and let the graph H be the union of the dense 
distance graphs of the pieces in V . Fakcharoenphol 
and Rao [FR06] devised an ingenious implementation 
of Dijkstra's algorithm | Dij59] that computes a shortest- 
path tree in H in time 0(\V(H)\ lg 2 n), where \V(H)\ 
is the number of nodes in H (i.e., the total number of 
boundary nodes in all the pieces in V). We will refer to 
this implementation as FR-Dijkstra. 

In fact, the proof of Fakcharoenphol and Rao's 
algorithm only relies on the property that the distances 
in each of the dense distance graphs given as input 
correspond to distances in a planar graph between a 
set of nodes that lie on a constant number of faces. It 
does not rely on any other properties of the r-division. 

FR-Dijkstra can be extended to the following set- 
ting (cf. BSWN10J). Let J be a planar graph. Let n' 
denote the number of nodes of J U H . We can compute 
shortest paths bxHUJ in 0(\H\ lg 2 \H\ + n' lg n 1 ) time. 
The edges of H are relaxed using FR-Dijkstra, while the 
edges of J are relaxed as in a traditional implementation 
of Dijkstra's algorithm using a heap. 

2.3.1 External Dense Distance Graphs. For a 

piece P, let G — P be the graph obtained from G by 
deleting the nodes in Vp — dP. The external dense 
distance graph of P, denoted by DDGg-p, is the 
complete graph on dP such that the length of an arc 



corresponds to the distance between its endpoints in 
G — P. External dense distance graphs were used 
recently in [BSWN10] . Computing the external dense 
distance graphs for all pieces in an r -division cannot be 
done efficiently using Klein's MSSP. The reason is that 
\V(G - P)\ can be 0(|V(G)|) for all pieces. Instead, 
the computation can be done in a top-down approach 
as follows (cf. BSWN10 ). Recall that an r -division 
is obtained as the set of pieces of the lowest level in 
a recursive division of the graph using Miller's simple- 
cycle separators. Consider the set Q of all pieces at all 
levels of the recursion (rather than just the set of pieces 
at the lowest level). Note that there are 0(lg|V(G)|) 
recursive levels since the size of the pieces decreases 
by a constant factor for every constant number of 
applications of Miller's cycle separator theorem. First, 
we compute DDGq for every piece Q £ Q. As explained 
above, this can be done in 0(\V{G)\lg\V{G)\) for all 
the pieces in a specific level of the recursion, for a total 
of 0{\V{G)\\g 2 \V{G)\) for all pieces in Q. Next, we 
consider a piece Q and denote the two subpieces of Q 
by Qi and Qi- DDGg-q 2 is obtained by computing 
distances in DDGq 1 U DDGq-q (see Figure [2]), using 
multiple applications of FR-Dijkstra, once for each node 
in 8Q 2 - This takes 0(\Q\ lg 2 \dQ\) time per piece, 
for an overall 0(\V(G)\ lg 2 |F(G)|) time for all pieces 
in a specific level. Since the number of levels is 
bounded by 0(lg | V(G)|), the entire computation takes 
0(|l/(G)|lg 3 |y(G)|)time. 




Figure 2: Pieces Q% and Qi in some level of the recursive 
application of Miller's cycle separator theorem. The 
piece Q is the union of Qi and Q2 ■ Distances in G — Q2 , 
the exterior of Q/2, are obtained by considering shortest 
paths in the interior of Q\ and in G — Q, the exterior 
of Q. 



3 A Linear-Space Distance Oracle. 

We first provide our linear-space data structure. The 
techniques employed here are reused in the subsequent 
sections, particularly in the cycle MSSP data structure 
(Section |4]). In the following we prove a more precise 
version of Theorem 11.31 



Theorem 11.31 For any directed planar graph G with 
non-negative arc lengths and for any constant e > 0, 
there is a data structure that supports exact distance 
queries in G with the following properties: the data 
structure can be created in time 0{n\gn), the space 
required is 0{n), and the query time is 0(n 1 / 2+e ). 

For non-constant e > 0, the preprocessing time is 
0(n lg(n) lg(l/e)), the space required is 0(nlg(l/e)) ; 
and the query time is 0(n 1 ^ 2+e + n 1//2 lg 2 (n) lg(l/e)). 

Our distance oracle is an extension of the oracle in 
Fakcharoenphol and Rao [FR06J. The main ingredients 
of our improved space vs. query time tradeoff are (i) 
using recursive r-divisions instead of cycle separators, 
and (ii) using an adaptive recursion]^] where the ratio 
between the boundary sizes of piece at consecutive levels 
is \pa. 

We split the proof into descriptions and analysis 
of the preprocessing and query algorithms. Let k = 

e(ig(i/e)). 

Preprocessing. We compute the recursive r- 
division of the graph with k recursive levels and values 
of {^i}f = o to be specified below. This takes 0{kn lgn) 
time. We then compute the dense distance graph for 
each piece. This is done for a piece P, with r nodes 
and 0(y /f r) boundary nodes on a constant number of 
holes, by applying Klein's MSSP algorithm [Kle05] as 
described in Section [2~3| Thus, all of the boundary-to- 
boundary distances in P are computed in 0(r lg r) time. 
Summing over all 0{n/ri) pieces at level i, the prepro- 
cessing time per level is 0(n\gri). The overall time to 
compute the dense distance graphs for all pieces over all 
recursive levels is therefore 0(knlgn). 

The space required to store DDGp is 0((^/r) 2 ) = 
O(r); summing over all pieces at level i we obtain space 
O(^Ti) = 0(n) per level; the total space requirement 
is O(kn). 

Query. Given a query for the distance between 
nodes u and v, we proceed as follows. For simplicity 
of the presentation, we initially assume that neither u 
nor v are boundary nodes. 

Let Pq be the level-0 piece that contains u. We 
compute distances from u in P - This is done in O(r ) 
time using the algorithm of Henzinger et al. [HKRS97 . 
Denote these distances by distp (u, w) for w € Pq- Let 
Hq denote the star graph with center u and leaves w € 
8Pq. The arcs of Hq are directed from u to the leaves, 
and their lengths are the corresponding distances in Po- 

Let S u be the set of pieces that contain u. Note 
that S u contains exactly one piece of each level. Let R u 

5 This appears to be the main difference to the distance oracle 
of Nussbaum |Nusll| . which uses a non-adaptive recursion. Using 
our adaptive recursion is crucial whenever S £ [n . lg lg n, n lg n). 



be the union of subpieces of every piece in S u . That 
is, R u = \Jp e s u {P' '■ P' is a subpiece of P}, Let H u 
be the union of the dense distance graphs of the pieces 
in R u . We use FR-Dijkstra (see Section [2~3| to compute 
distances from u in H u U Ho. Observe that any shortest 
path from u to a node of H u can be decomposed into 
a shortest path in P$ from u to dPo and shortest paths 
each of which is between boundary nodes of some piece 
in R u . Since all u-to-dPo shortest paths in Pq are 
represented in H , and since all shortest paths between 
boundary nodes of pieces in R u are represented in H u , 
this observation implies that distances from u to nodes 
of H u in H u U H are equal to distances from u to nodes 
of H u in G. We denote these distances by dist^^it, w) 
for nodes w € H u . 

We repeat a similar procedure for v (reversing the 
direction of arcs) to compute distG(w, v), the distances 
in G from every node w € H v to v. 

Let P uv be the lowest-level piece that contains both 
u and v. Assume first that P uv is not a level-0 piece. 
Let P u (P v ) be the subpiece of P uv that contains u (v). 
Since P uv is both in S u and in S v , both P u and P v are in 
R u as well as in R v . This implies that we have already 
computed distc(u,?iO and dist^io, i>) for all w £ dP u . 
Since we have assumed that P uv is not a level-0 piece, 
the shortest u-to-v path must contain some node of 
dP u . Therefore, the u-to-v distance can be found by 
computing 

min distc(u, w) + distc(w, v). 

If P uv is a level-0 piece, then P uv — P , and the 
u-to-v distance can be found by computing 



min < distp (u, w), min {dist^n, w) + distc(w, v)} 



The case when u or v are boundary nodes is 
a degenerate case that can be solved by the above 
algorithm. Let Q u be the highest-level piece of which u 
is a boundary node. We have the preprocessed distances 
in Q u from u to all other nodes of dQ u . Therefore, it 
suffices to replace S u above with the set of pieces that 
contain Q u as a subgraph in order to assure that H u 
is small enough and that the distances computed by 
the fast implementation of Dijkstra's algorithm are the 
distances from u to nodes of H u in G. 

Query Time. Computing the distances distp (-, •) 
takes O(ro) time. Let |V(ff u )| denote the number 
of nodes of H u . The FR-Dijkstra implementation 
runs m 0(\V{H u )\lg 2 n) time. It therefore remains 
to bound |V(ffu)|. Let Pi be the level-i piece in S u . 
Pi has 0(^7) subpieces, each with 0(y/ri_i) boundary 
nodes. Therefore, the contribution of Pi to |y(i? u )| is 



Q{ ^* i ). The total running time is therefore 

Recall that = |V^(G)| = n, and set r$ = y/n. For 
i = 1 . . . k — 1 we recursively define n so as to satisfy 



— vn. 



This implies 



iii i_i 
r\ = rt 2 4 = n 4 

1+3 jl 



n 2 ^ is = n is 



r fc _! = n 2k . 
The total running time is thus bounded by 



O V^ + lg 2 ™ (k-l)Vn + 



V n 2k / 
O ( (ky/n + n^ + i^TT j i g 2 n 



< 



(3.1) 



By setting k — 0(lg(l/e)) we obtain the claimed 
running times. 

4 A Cycle MSSP Data Structure for Planar 
Graphs. 

In this section we provide our main technical tool. We 
prove Theorem 1 1.4 1 which we restate here. 



Theorem 11.41 Given a directed planar graph G on n 
nodes and a simple cycle C with c = 0{y/n) nodes, there 
is an algorithm that preprocesses G in 0(nlg 3 n) time 
to produce a data structure of size O(nlglgc) that can 
answer the following queries in 0{c lg 2 clglgc) time: for 
a query node u, output the distance from u to all the 
nodes of C . 

Comparison with Klein's MSSP data struc- 
ture Our data structure can be seen as an alternative 



to Klein's MSSP data structure (see Section 2.2) with 
two main advantages (which we exploit in Section [5]) : 

• our data structure can handle queries to any not- 
too-long cycle as opposed to a single face (the 
crucial difference and difficulty is that shortest 
paths may cross a cycle but not a face), 

• the space requirements are only O(nlglgc) (we in- 
ternally rely on the data structure in Theorem [L3j 
so even 0(n) is possible at the cost of increasing 
the query time) as opposed to O(nlgn), 



and three main disadvantages: 

• our data structure cannot efficiently answer queries 
from u to a single node on the cycle G; such a 
query requires the same time as computing the 
distances from u to all the nodes on G (in many 
existing applications, this disadvantage is not really 
problematic, since MSSP is used for all nodes on 
the face anyway), 

• our data structure requires amortized time 
0(lg 2 clglgc) per node on the cycle (as opposed 
to O(lgn)), which is slower for long cycles, and 

• the preprocessing time of our data structure is 
0(nlg 3 n) as opposed to 0(n lgn). 

Let G be an embedded planar graph. 

Preprocessing. Let Go be the exterior of G. That 
is, the graph obtained from G by deleting all nodes 
strictly enclosed by C . Consider C as the infinite face 
of Go- Similarly, let Gi be the interior of G. Namely, 
the graph obtained from G by deleting all nodes not 
enclosed by C. Consider G as the infinite face of G\. 
Note that we have reduced the problem from query 
nodes on a cycle G to query nodes on a face. The query 
algorithm is supposed to handle paths that cross G. The 
preprocessing step consists of the following: 

1. Computing DDGc and DDGq-c- This can be 
done in 0((n + c 2 )lgn) time using Klein's MSSP 
algorithm [Kle05j . Storing DDG C and DDGg-c 
requires 0(c 2 ) = 0(n) space. 

2. Computing an r-division of Gj (i € {0,1}) with 
r = c 2 . Each piece has 0(c 2 ) nodes and 0(c) 
boundary nodes incident to a constant number of 
holes. Consider all the nodes of G as boundary 
nodes of every piece in the division. Note that each 
piece still has 0(c) boundary nodes. This step takes 
O(nlgn) time. 

3. Computing, for each piece P, a recursive r-division 
of P as the one in the preprocessing step of the 
oracle in Section [3] (Theorem 1.3) with e = 1/lgc. 
That is, the number of levels in this recursive r— 
division is k — O(lglgc). The top-level (level-fc) 
piece in this recursive division is the entire piece 
P. In the description in Section [3j the top-level 
piece is the entire graph and therefore it has no 
boundary nodes. Here, in contrast, we consider 
the boundary nodes of P as boundary nodes of the 
top-level piece in the decomposition (and thus, as 
boundary nodes of any lower-lever piece in which 
they appear) . This does not asymptotically change 
the total number of boundary nodes at any level 



of the recursive decomposition since P has 0(c) 
boundary nodes, and every level of the recursive 
decomposition consists of a total of 0(c) boundary 
nodes. The time to compute the recursive r- 
division for all pieces is bounded by 0(?ilg 2 n). 

4. Computing, for each piece P, the dense distance 
graph for each of the pieces in the recursive de- 
composition of P. Let Hp denote the union of 
the dense distance graphs for all the pieces in 
the recursive decomposition of P. As discussed 
in Section |3j the space required to store Hp 
is 0(|P|lge) = 0(|P|lglgc). Using the meth- 
ods presented in Section [31 computing Hp takes 
0(|P|lg|P|lglgc) = O(cMgclglgc). Thus, the 
total space and time required over all pieces P is 
O(nlglgc) and 0(n lgclglgc), respectively. 

5. Computing, for each piece P, the dense distance 
graph DDGd-p ■ Recall that we consider the 
nodes of G as boundary nodes of every piece in 
the division. These dense distance graphs can 
be computed as described in Section |2.3.1| As 
shown there, the entire computation (for all pieces 
combined) takes 0(nlg 3 n) time. 

The time required for the preprocessing step is therefore 
0(nlg 3 n) and the space required is 0(n lglgc). 

Query. When queried with a node u, the data 
structure outputs the distances from u to all the nodes 
of G. We describe the case where u is not enclosed by G. 
In this case, we use the dense distance graphs computed 
in the preprocessing step for Go- The symmetric case 
is handled similarly, by using the dense distance graphs 
computed for G\. 

Let P be the piece in the r-division of Go to which 
u belongs. Recall that P consists of 0(c 2 ) nodes. 
Consider the recursive J — division of P computed in 
item [3] of the preprocessing stage. Let P be the level-0 
piece of P that contains u. Pq consists of 0( vc 2 ) = 0(c) 
nodes0 

We first compute }HKRS97j . in 0(c) time, the 
distances from u to all nodes of Po, and store them 
in a table distp . We then compute, using FR-Dijkstra, 
the distances from u in the union of the following dense 
distance graphs (see Figure [3]): 

1. H , the star graph with center u and leaves <9Po- 
The arcs of H are directed from u to the leaves 
and their lengths are the corresponding lengths in 
dist p . 

"Here, as in Section [3] we assume that u is not a boundary 
node in the recursive i — division of P. The case where u is such a 
boundary node is degenerate, see Section [3] 



Figure 3: A schematic diagram showing the various subgraphs whose dense distance graphs are used in a query 
to the cycle MSSP data structure. The cycle C is double-lined. The interior of C is the subgraph G\. The query 
node u is indicated by a small solid circle. The piece P in the r-division of the exterior of C (Gq) is shown as a 
grey region with solid boundary. The boundaries of the pieces whose dense distance graphs are in H u are shown 
as dotted lines (one level) and dashed-dotted lines (another level). Po is the smallest piece of P that contains u. 
Any shortest path from u to C can be decomposed into a shortest path from u to 8Pq followed by shortest paths 
between nodes on the boundaries shown in the figure. 



2. H u , the subset of dense distance graphs in Hp that 
correspond to pieces in the recursive decomposition 
of P that contain u and their subpieces. These 
dense distance graphs are available in H u . 

3. DDGo -p 

4. DDG C 

Note that the first two graphs are the analogs of H and 
H u from Section [3] 

Distances from u in the union of the above graphs 
are equal to the distances from u in G. This is true since 
any u-to-C shortest path can be decomposed into (i) a 
shortest path in P from u to 8Pq, (ii) shortest paths 
each of which is a shortest path in Q between boundary 
nodes of Q for some piece Q in the recursive r-division 
of P that is represented in H u , (Hi) shortest paths in 
Gq — P between nodes of dPUC, and (iv) shortest paths 
in the interior of C between nodes of C . 

To bound the running time of FR-Dijkstra we need 
to bound the number of nodes in all dense distance 
graphs used in the FR-Dijkstra computation. Hq has 
0(yfc) nodes. The analysis in Sectio n |3| shows that the 
graphs in the set H u consist of Q(^/ |P| lg lg \ P\) nodes 
(substitute k = l/lg|P| in eq. jpU])). DDG Go ^ p has 
0(c + y/\P\) = 0(c) nodes, and DDGc has c nodes. 
Combined, the running time of the invocation of FR- 



Dijkstra is bounded by 0(c lg 2 clglgc). This dominates 
the O(c) time required for the computation of distp , so 
the overall query time is 0{c lg 2 clglgc), as claimed. 

5 Distance Oracles with Space S G [nlglgn,n 2 ]. 

In this section we prove Theorem |1.1| Using our 
new cycle MSSP data structure, the proof is rather 
straightforward. 

Theorem 11.11 Let G be a directed planar graph on 
n vertices. For any value S in the range S € 
[nlglgn, n 2 ], there is a data structure with preprocess- 
ing time 0(S lg 3 n/ lglgn) and space O(S) that answers 
distance queries in 0{nS-^ 2 \g 2 n\g 3/2 \gn) time per 
query. 

Proof. Let r := (n 2 lg\gn)/S. Note that r e [lglgn, n] 
for any S G [nlglgn, n 2 ]. 

Preprocessing. We start by computing an r— 
division. Each piece has 0(r) nodes and 0(y/r) bound- 
ary nodes incident to a constant number of holes. For 
each piece P we compute the following: 

1. A distance oracle as in Theorem 11.31 with e = 
1/lgr. This takes 0(r lgrlglgr) time and 
O(rlglgr) space. 



2. For each hole of P (bounded by a cycle in G) 
we compute our new cycle MSSP data structure^ 
Since the number of holes per piece is constant, this 
requires requires 0(nlg 3 n) time and O(nlglgn) 
space per piece. 

Summing over all pieces, the preprocessing time is 
0(Slg 3 n/lg\gn) and the space needed is O(S). 

Query. Given a pair of nodes (s,t), we compute 
a shortest s-t path as follows. Assume first that s 
and t are in different pieces. Let P denote the piece 
that contains s and let dP denote its boundary. We 
compute the distances in G from dP to t using the cy- 
cle MSSP data structures. These distances can be ob- 
tained in time 0(\dP\ lg 2 nlglgn) = 0(^/r\g 2 nlglgra). 
Analogously, we compute the distances in G from s 
to dP. It remains to find the node p £ dP that 
minimizes da{s,p) + dc(p,t), which can be done in 
0{\dP\) = 0(y/r) time using a simple sequential search. 

If s and t lie in the same piece, the shortest s-t 
path might not visit dP. To account for this, we query 
the distance oracle for P, which takes 0(y/r lg 2 rlglgr) 
time. We return the minimum distance found. □ 

k— many distances As a consequence, we also 
obtain an improved algorithm for fc-many distances, 
whenever k — f2(-y/n/lglgn). 



Proof, [of Theorem 1.2 For some value of r to be spec- 
ified below, we preprocess G in time 0((n 2 /r) lg 3 n), 
and then we answer each of the k queries in time 
0(^/r\g 2 rlglgr). The total time is 0((n 2 /r) lg 3 n + 
ky/r lg 2 r lg lg r). This is minimized by setting r = 
n 4 / 3 fc~ 2 / 3 (lgn/lglgn) 2 / 3 . Note that r = 0(n) since 
k = Q(y/n\gn/ \g\gn). The total running time is thus 
0((fcn) 2 / 3 (lgn) 7 / 3 (lglgn) 2 / 3 ). □ 

Comparison and Discussion. The query time 
of our data structure is at most 0(y/r lg 2 rlglgr), 
which, in terms of S, is OfjiS^ 1 / 2 lg 2 n lg 3 ^ 2 lgn). 
Let us contrast this with Cabello's data struc- 
ture |Cab06j that, for any S £ [n 4 / 3 lg 1/3 n, n 2 } has 
preprocessing time and space O(S) and query time 
0(nS'~ 1 / 2 lg 3 / 2 n). In our construction, we sacrifice 
a factor of 0(y / \gn(\g lgn) 3 ) in the query time but 
we gain a much larger regime for S. For the range 
S € [w(rtlgn/lglgn),o(n 4/ ' 3 lg 1 ^ 3 n)], only data struc- 
tures of size O(S) with query time 0(n 2 / S) had been 
known Dji96] (see also Figure ll}. 

To conclude this section, let us observe what hap- 
pens when we gradually decrease the space requirements 
S from n 2 down to n, and, en passant, let us pose some 



open questions. For quadratic space (or even slightly 
below [WNlOaj), we can obtain constant query time. 
As soon as we require the space to be 0(n 2 ~ e ) for some 
e > 0, the query time increases from constant or poly- 
logarithmic to a polynomial. It is currently not known 
whether 0(1) is possible or not — the known lower 
bounds |SVY09[ IPR10] on the space of distance oracles 
work for non-planar graphs only. Further restricting the 
space, as long as the space available is at least il(nlgn), 
the data structure can internally store MSSP data struc- 
tures. The query time for this regime of S can actually 
be made slightly faster than what we claim in Theo- 
rcm |l.l| (by avoiding the 0(lglgn)-factor due to the re- 
cursion needed in Theorem 11.31 and its manifestation in 
the cycle MSSP data structure). The data structure of 
Nussbaum [Nusllj obtains a "clean" tradeoff with query 
time proportional to 0(n/y/S) for S > n 4 / 3 (without 
logarithmic factors in the query time, currently at the 
cost of a slower preprocessing algorithm). The obvious 
open question is whether it is possible to obtain a data 
structure with space O(S) and query time 0(n/\/S) for 
the whole range of S £ [n, n 2 ] and without substantial 
sacrifices with respect to the preprocessing time. An- 
other open question is whether it is possible to improve 
upon this tradeoff. Note that, for quadratic space, an 
improvement of an (almost) logarithmic factor is possi- 
ble [WNlOaj . 

As soon as we require the space to be o(n lgn), 
we cannot afford to store the MSSP data structure 
anymore and we are currently forced to rely on the cycle 



MSSP data structure (Theorem 1.4). The query time 



Y Note that for S = o(n lg n) we cannot even afford to store 
Klein's MSSP data structure |Kle05| . 



increases to the time bound claimed in Theorem 11.11 
When we further restrict the space to o(nlglgn), say 

5 = 0(nlg(l/e)) for some e > 0, the query time 
increases to 0(n 1//2+e ). A different recursion (maybe 
a la [HKRS97 ) could potentially reduce this to 0(y/n). 

Let us briefly consider the additional space of the 
data structure, assuming that storing the graph is free. 
It is known that, for approximate distances, a query 
algorithm can run efficiently using a data structure that 
occupies sublinear additional space KKSl 1|. For exact 
distances and sublinear space, nothing better than the 
linear-time SSSP algorithm HKRS97 is known. 

6 Distance Oracles with Query Time 
Quasi-Proportional to the Shortest-Path 
Length 

We use our new cycle MSSP data structure to prove 
Theorem |1.5| which states that there is a distance 
oracle with query time proportional to the shortest-path 
length. We actually prove two versions, the stronger one 
being a distance oracle with query time proportional 
to the minimum number of edges (hops) on a shortest 



path. For the stronger version, we need the following 
assumption, which essentially means that approximate 
shortest paths do not use significantly fewer edges. 

Assumption 6.1. Let hc(s,t) denote the number of 
edges on a minimum-hop shortest-path. For some 
constant e > 0, all s-t paths of length at most (1 + 
e)d(}(s, t) have Cl(JiG(s,t)) edges. 

We restate Theorem 1 1 . 5 1 and its stronger variant. 



Theorem 11.51 For any planar graph G with edge 
lengths > 1 there is an exact distance ora- 
cle using space 0(n\gn\glgn) with query time 
0(min {£ lg 2 £ lg \g£, \/n\g 2 n\) for any pair of nodes 
at distance £. The preprocessing time is bounded 
by 0(n 1+e ) for any constant e > 0. 

Furthermore, if Assumption \ 6. 1\ holds for G, the 
query time is at most 0(min{Mg h lg lg h, -^/nlg 2 n}) 
for any pair of nodes (u, v) at hop-distance h — 
h G (u,v). 

Our main ingredient is a distance oracle for planar 
graphs with tree-widtl^w. 

Theorem 6.1. Let G be a planar graph on n vertices 
with tree-width w. For any constant e > and for 
any value S in the range S € [nlglgw,n 2 ], there is 
a data structure with preprocessing time 0(S\g 2 n + 
n 1+e ) and space O(S) that answers distance queries in 
0(min{nS'~ 1 / 2 lg 2 ' 5 n, w lg 2 w lg lg w}) time per query. 

Note that for 5* superlinear in n but less than roughly 
nw, the oracle cannot make any use of the additional 
space available and the query algorithm runs in time 
proportional to w (up to logarithmic factors). Any 
application should use either space close to linear or 
more than fl(nw). 

Using Theorem |6.1| the proof of Theorem |1.5| 
boils down to a combination of slicing (see for exam- 
ple |Bak94[lKle08] y local tree-width |Epp00[ IDH04] . and 
scaling. 



Proof, [of Theorem 1.5 We repeatedly use Theorem 6.1 
for different subgraphs as follows. 

Slicing and local tree-width. We use standard 
techniques |Bak941 IKle08j and research on the related 
linear local tree-width property EppOO , DH04] to slice a 
planar graph into subgraphs of a certain tree-width: 

• We compute a breadth-first search tree in the 
planar dual rooted at an arbitrary face. 



• The subgraph induced by the nodes at depth d £ 
[d 1 ,d 2 ] has tree- width 0{d 2 - d x ) |Epp00[ [DHOU . 

If we were only interested in paths using at most w 
edges, we could (i) cut the graph into slices of depth w, 
and (ii) check the union of any two consecutive levels 
containing both endpoints. It is straightforward to see 
that each path on w edges lies completely within one of 
these unions. 

Scaling. We apply the slicing step described in the 
previous paragraph for different scales. 

• For every integer i > with 2* < y/n, we slice the 
graph into subgraphs GJ at depth r = 2 4 ; here, G] 
denotes the graph induced by the nodes adjacent 
to all the faces at depth in [jr, (j + l)r). Every 
subgraph G* has tree- width at most 0(r). 

• For any two consecutive G* , G* +1 we compute a dis- 
tance oracle of size 0(\V(G\ U G) +1 )\ lg lg \V(G) U 
Gj +1 )\) with query time 0(r lg 2 rig lgr) as in The- 



orem |6T] Since each node is in at most two graphs 
Gj and since each G* participates in at most two 
distance oracles, the total size of all these distance 
oracles per level i is 0{n lglgn). The total size of 
our data structure is thus O(nlgnlglgn). 

Which Scale? Let the smallest number of edges 
(or hops) on any shortest path from u to v be h = 
hc(u, v). At query time we use an approximate distance 
oracle to determine the right scale. In the preprocess- 
ing algorithm, we also precompute the approximate dis- 
tance oracle of Thorup |Tho04| for e = 1/2. This oracle 
can be computed in 0(ne _1 lg 3 n) time, it uses space 
0(ne _1 lgn), and it answers (1 + e)-approximate dis- 



s See |Hal76l IRS86I IAP89| for definitions and much more on 
tree- width. 



tance queries in time 0(l/e). If Assumption 6.1 holds, 
we instead use that value of e in the construction of 
Thorup's distance oracle. The asymptotic space con- 
sumption is not increased. 

Query Algorithm. At query time, given a pair of 
nodes (it, v) at distance £, we need to find a level that 
contains a shortest path. We query the approximate 
distance oracle in time 0(l/e) to obtain an estimate 
for £. Let £ denote this estimate. We then execute one 
of the following search algorithms. 

In the case that Assumption |6.1| does not hold, we 
directly query level i for the smallest i with 2 l > £. Since 
all the edge weights are at least 1, any path of length £ 
has at most £ edges and is thus contained in some graph 
Gj at level i. Since graphs at level i have tree-width 
0(2^, and since 2 l = 0(1) = 0(1), the running time of 
the query algorithm is 0(£\g 2 £\g\g£) as claimed. 

If Assumption |6.1| does hold, we search the data 
structure level by level with increasing i until the first 



time a distance at most £ is found. By Assumption 6.1 
we know that any u-to-v path of length < (1 + t)t uses 
at least c • h G (u,v) edges for some constant c < 1. We 
therefore search the next 1 — lg 2 c levels to ensure that 
we find a shortest path. The running time is a geometric 
sum, which is dominated by the time to search the last 
level with tree-width 0(hc(u, v)). □ 

6.1 Distance Oracles for Planar Graphs with 
Tree- width o(y / n). In this section we prove Theo- 
rem 6.1 There exist efficient distance oracles for (not 
necessarily planar) graphs with tree- width w. Chaud- 
huri and Zaroliagis [CZOOJ provide a distance oracle 
that uses space 0(w 3 n) and answers distance queries in 
time 0(u> 3 a(n)), where a(n) denotes the inverse Ack- 
crmann function. Farzan and Kamali |FK11) provide 
a distance oracle that uses space 0(wn) and answers 
distance queries in time 0(w 2 lg 3 w). 

In our application, the tree-width w may be non- 
constant up to O(yfn). For this reason we cannot use 
their distance oracles. In the following we improve upon 
their results for the special case where graphs are further 
assumed to be planar: we can obtain space 0(n lglgiw) 
and query time 0(w\g 2 wlglg w) using our distance 
oracle. 

In some sense our proof can be seen as an improve- 
ment over Djidjev's data structure with space O(S) 
and query time 0{n 2 /S). His data structure works on 
any graph with recursive balanced separators of size 
0(y/n) and, furthermore, similar results can be ob- 
tained for graphs with separators of general size f(n) — 
o(n) [Dji96[ Sections 3 and 4]. The data structure 
with the better trade-off (space O(S) and query time 
0{n/\/~S)) however only works for planar graphs |Dji96| 
Section 5], since it exploits the Jordan Curve Theo- 
rem. We can now exploit the Jordan Curve Theorem 
for planar graphs with smaller separators by observing 
that, for planar graphs with smaller tree- width (o(y / ri)), 
the size of the Jordan Curve separating the inside from 
the outside decreases proportionally |ST941 DPBF10 . 
These separators are referred to as sphere-cut separa- 
tors |ST941 IDPBFlOj and they can be found efficiently. 



Lemma 6.1. (Gu and Tamaki |GT09i Theorem 2]) 
For any constant e > and for any biconnected vertex- 
weighted planar graph on n nodes with tree-width w, 
there is an 0{n l+,i )-time algorithm that finds a non- 
self-crossing cycle C of length 0(w/e) such that any 
connected component of G\C has weight at most 3/4 
the total weight. 

Proof. The algorithm computes a branch decomposition 
as in Gu and Tamaki [GT091 proof of Theorem 2]. It 
is known that branch-width and tree-width are within 



constant factors |RS91j . 

A branch decomposition RS91 of a graph G is a 
ternary tree T whose leaves correspond to the edges 
of G. Removing any tree edge e <E T creates two 
connected components Ti,T2 C T, each of which 
corresponds to a subgraph G,; of G, induced by the 
edges corresponding to the set of leaves of T in Tj . In 
the branch decomposition tree, we may thus choose the 
edge e* that achieves the best balance among all the 
edges e e T. 

In the algorithm of Gu and Tamaki, each edge of the 
branch decomposition tree corresponds to a self-non- 
crossing closed curve passing through 0(w/e) nodes of 
G (as in sphere-cut decompositions) that encloses Gi 
and does not enclose Gi (or vice versa). 

Since the degree of each node is at most three, it is 
possible to choose an edge e* such that the total weight 
strictly enclosed by the non-self-crossing cycle G e * and 
the total weight strictly not enclosed by G e » are each at 
most a 3/4-fraction of the total weight. The cycle G e * 
is our balanced separator of length 0(w/e). □ 

For an optimal sphere-cut decomposition, the 
fastest algorithm runs in cubic time [ST94 , GT08 . For 



our purposes, the constant approximation in Lemma 6.1 
suffices. 

We define a variant of the r-division as in Sec- 
tion |2.1[ wherein we use either Miller's cycle separator or 
sphere-cut separators, whichever is smaller. An [r, w]- 
division of G with tree-width w is a decomposition into 
<d(n/r) edge-disjoint pieces, each with 0(r) nodes and 
0(m\rt{y/r, w}) boundary nodes. 

Lemma 6.2. For any constant e > and for any planar 
graph on n nodes with tree-width w, an [r, u>]-division 
can be found in time 0(n 1+e lg n + n lg r + nr^ 1 / 2 lg n) . 

Proof. Fix e to be the desired constant > 0. The 
procedure for obtaining the [r, u>]-division is the one 
described in [WNlObj . but using either Miller's cycle 
separator or sphere-cut separators, whichever is smaller. 
If yfr = 0{w) then since the separators we use are at 
least as small as the separators assumed in the proof 
of |WN10bl Lemmata 2 and 3], the lemmas apply and 
the resulting decomposition has 0(n/r) pieces, each 
with 0(r) nodes and 0(^/r) boundary nodes. 

If y/r = w(w) then first apply the separator the- 
orem as in the procedure for obtaining a weak r— 
division WNlObl Lemma 2] until every piece has size 
0{r). Note that since \/r — w(w), we always use sphere- 
cut separators, whose size is 0(w) (hiding a linear factor 
of 1/e), regardless of the size of the piece we separate. 
Therefore, by [WNlOb, Lemma 2], the number of pieces 
is 0(n/r), and the total number of boundary nodes is 



0(n/y/r). However, here we need a tighter bound on 
the total number of boundary nodes. In the following 
we show that the total number of boundary nodes is 
0{nw/r). 

Consider the binary tree whose root corresponds 
to G, and whose leaves correspond to the pieces of 
the weak r-division. Each internal node in the tree 
corresponds to a piece in the recursive decomposition 
that either consists of too many nodes or contains 
too many holes and is therefore separated into two 
subpieces using a sphere-cut separator. As in the proof 
of [WNlOb, Lemma 2], for any boundary vertex v in 
the weak r-division, let b(v) denote one less than the 
number of pieces containing v as a boundary vertex. 
Let B{n) be the sum of b(v) over all such v. Every time 
a separator is used, at most cw boundary nodes are 
introduced for some c > 0, and the piece to which these 
nodes belong to is split into two pieces, so the number of 
pieces to which these nodes belong to increases by one. 
Therefore, B(n) is bounded by the number of internal 
nodes times cw. Since the tree is binary, the number of 
internal nodes is bounded by the number of leaves, so 
B(n) = 0{nw/r), which shows that the total number 
of boundary nodes in this weak r-division is 0(nw/r). 

Next, consider the procedure in [WNlOb, Lemma 3], 
which further divides the pieces of the weak r-division 
to make sure that the number of boundary nodes in each 
piece is small enough. In our case we want to limit the 
number of boundary nodes per piece to at most c'w for 
some d (which we set > c below). Let ti denote the 
number of pieces in the weak r-division with exactly i 
boundary nodes. Note that Yli^i = J2vev B (^( v ) + 1)' 
where Vb is the set of boundary vertices over all pieces 
in the weak r-division. Hence, XL^tj < 2B(n), so by 
the bound on B{n), J^viij — 0(nw/r). 

In the weak r-division, consider a piece P with i > 
c'w boundary vertices. When the above procedure splits 
P, each of the subpieces contains at most a constant 
fraction of the boundary vertices of P. Hence, after 
di/(c'w) splits of P for some constant d, all subpieces 
contain at most c'w boundary vertices. This results 
in at most 1 + di/(dw) subpieces and at most cw new 
boundary vertices per split. We may choose d to be 
sufficiently larger than c. The total number of new 
boundary vertices introduced by the above procedure 
is thus 

2_^cw(di/(dw))ti < itj = 0(nw/r) 

i i 

and the number of new pieces is at most 
Y^{di/{cw))U = 0(n/r). 



Hence, the procedure generates an [r, w]-division. Since 
a sphere-cut separator can be found in 0(n 1+e ) time, 
a weak r-division can be found in 0(n 1+e \gn) — 



Ob 



time, and refining it into an [r, w] -division can 



be done within this time bound. 



□ 



Proof, [of Theorem 6.1 ] The data structure for planar 



graphs with smaller tree-width is essentially the same 
as the data structure for general planar graphs as 
described in Section [5] The main difference is that, 
when computing a cycle separator, instead of Miller's 
algorithm to find a cycle separator of length 0(y/n), we 
use sphere-cut separators and the corresponding [r, w\- 
division as in Lemma 16.21 □ 

7 Conclusion. 

We introduce a new data structure to answer distance 
queries between any node v and all the nodes on a 
not-too-long cycle of a planar graph. Using this tool, 
we significantly improve the worst-case query times for 
distance oracles with low space requirements S (down 
from 0(n 2 /S) to 0{yjn 2 / S)). Furthermore, for linear 
space, we improve the query time down from 0(n) 
to 0(n 1 / 2+e ) using adaptive recursion. We also give 
the first distance oracle that actually exploits the tree- 
width of a planar graph, particularly if it is o{^/n) 
and, as an application, we give a distance oracle whose 
query time is roughly proportional to the shortest-path 
length. Similar behavior of practical methods had been 
observed experimentally before but could not be proven 
until the current work. 

An interesting and important open question is 
whether there is another tradeoff curve below the space 
O(S) and query time 0(n/V~S) curve. 
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